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Abstract 

We explore the compact U(l) lattice gauge theory with staggered fermions and 
gauge field action — X)p[/3cos(0p) +7cos(20p)], both for dynamical fermions and 
in the quenched approximation. (0p denotes the plaquette angle.) In simulations 
with dynamical fermions at various 7 < —0.2 on 6^ lattices we find the energy gap at 
the phase transition of a size comparable to the pure gauge theory for 7 < on the 
same lattice, diminishing with decreasing 7. This suggests a second order transition 
in the thermodynamic limit of the theory with fermions for 7 below some finite 
negative value. Studying the theory on large lattices at 7 = —0.2 in the quenched 
approximation by means of the equation of state we find non-Gaussian values of the 
critical exponents associated with the chiral condensate, /? ~ 0.32 and 5 ~ 1.8, and 
determine the scaling function. Furthermore, we evaluate the meson spectrum and 
study the PCAC relation. 



1 Introduction 

Many years of nonperturbative lattice studies have not yet clarified the prop- 
erties of the strongly coupled QED. Concerning the noncompact lattice formu- 
lation with staggered fermions the properties of the 2'^'^ order phase transition 
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between the phase with broken chiral symmetry at strong couphng and the 
symmetric weak couphng (Coulomb) phase are stih controversial (see [1-3] 
and references therein). So is the existence of the conjectured non-Gaussian 
fixed point which might allow to avoid triviality of QED in the limit of in- 
finite cutoff. The compact strongly coupled QED (U(l) gauge theory) with 
staggered fermions was found to have a phase with broken chiral symmetry, 
too [4-7]. But the corresponding phase transition to the symmetric phase has 
not been considered of interest for the continuum theory, as this transition was 
assumed to be of 1®* order because of an observed two-state signal on finite 
lattices [8,9]. 

There are two reasons to reconsider the order of the phase transition in com- 
pact QED. Firstly, a suitable modification of the lattice action might change 
the order. Such a phenomenon occurs in the pure compact U(l) lattice gauge 
theory with the extended Wilson action 

'5t/ = -E[/3cos(ep)+7cos(2ep)]. (1) 
p 

(0p is the plaquette angle proportional to F^^. Further details are given in 
sec. 2.) There the analysis of the decreasing two-state signal [10,11] suggested, 
that the order of the transition line (3 = Pdl) changes from the I''* to the 2"*^ 
with decreasing 7. In the full theory with fermions of bare mass mo (in lattice 
units) the two-state signal was found to decrease with decreasing 7 for fixed 
lattice size and mo, too [9,12]. 

Secondly, the recent investigation of the pure U(l) gauge theory with the ex- 
tended Wilson action [13-15] and the Villain action [16] on finite lattices with 
sphere-like topology revealed that the two-state signal might in part be due 
to finite size effects of topological origin. In some range of couplings, when the 
two-state signal is rather weak, it vanishes when the lattices with sphere-like 
topology instead of the usual toroidal ones arc used, strongly suggesting 2"*^ 
order in the thermodynamic limit. The two-state signal may be related to cer- 
tain gauge field configurations, e.g. monopole loops winding around a toroidal 
lattice. Since the simulations of the full compact QED were performed on cu- 
bic lattices with periodic boundary conditions for gauge fields, the same phe- 
nomenon should be expected also here, quite independently of the fermions. 
Thus it may well be that in the thermodynamic limit the two-state signal 
again vanishes and the phase transition is actually of 2°^^ order. 

In this paper we therefore hypothesize that in compact QED with the extended 
Wilson action and light staggered fermions a part of the phase transition 
surface (parametrized by 7 and mo) is of 2"^^ order. It would be challenging 
to check this conjecture in simulations with dynamical fermions on lattices 
with sphere-like topology. Unfortunately, putting fermions on such lattices is 
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demanding, requiring various preparatory studies, e.g. with free fermions. So 
presently we are not in the position to make such a direct test. 

On the usual toroidal lattices of fixed sizes one can at least investigate, how 
the two-state signal depends on 7 and mo and compare its size to that of the 
pure U(l) gauge theory. Here we report on such a systematic study of the two- 
state signal on 6'^ lattice. The size of the two-state signal can be conveniently 
characterized by the difference ("gap") Aep between the mean values ep of 
cos(6p) in each of the coexistent states. 

A study of the 7-dependence of Aep at fixed mo was made by Okawa [12] 
some time ago. He suggested that at 7 = —1.3 and mo = 0.1 the two-state 
signal is already absent and the transition is therefore of 2^^^^ order there. Our 
present study with larger statistics and a modern hybrid Monte Carlo fermion 
algorithm shows that the two-state signal persists on the toroidal 6*^ lattice 
even at 7 = —1.5 (simulations at still smaller 7 are prohibitively expensive). 
However, Aep is smaller than in the pure U(l) theory at 7 ~ —0.5, where 
recent work on spherical lattices [14,15], as well as earlier investigations on 
toroidal lattices [11] strongly suggest 2^^ order. 

At fixed 7 the gap Aep increases as mo decreases. Since in the continuum 
limit fermions of finite mass can be obtained only when mo 0, the value 
of Aep at mo = must be considered. At 7 = —0.2 we find that Aep, 
extrapolated to this limit, is of the same magnitude as in the pure gauge 
theory at 7 = 0. Numerical evidence on spherical lattices suggests 2°^^ order in 
the latter case [13-15]. Thus our results of dynamical fermion calculations on 
6'^ lattice, described in sec. 3, support the hypothesis that at some negative 7 
the transition is indeed of 2^^^^ order even for light fermions. 

If so, then it is justified to consider the continuum limit and investigate its 
physical properties. We do first steps in this direction, using the quenched 
approximation. As described in sec. 4, we have studied on toroidal lattices 
at 7 = —0.2 the chiral condensate, susceptibilities, equation of state and 
the meson masses. The data have been obtained in Monte Carlo runs for a 
simultaneous investigation of the gauge-ball spectrum in the pure U(l) gauge 
theory in both phases [17]. For that purpose we have accumulated a large 
number of gauge field configurations at 17 /3-points on lattices up to 20^40 
(further details are given in [17]). All these points are close to the phase 
transition, but outside the tunneling region where the two-state signal occurs. 
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2 Action and phase diagram 



Compact QED with extended Wilson action for gauge field and staggered 
fermions on the lattice is defined by the action 



S^S^ + Su, (2) 

Sx-T.X.M,,yXy (3) 

1 ^ 

'5c/ = -E[/^cos(ep)+7cos(2ep)]. (5) 



Here, 0p G [0,27r) is the plaqucttc angle, i.e. the argument of the product 
of U(l) link variables around a plaquette P, and P and 7 are the single and 
the double charge representation couplings, respectively. Taking Qp = o^gF^^, 
where a is the lattice spacing, and (3+4'-/ = l/g"^, one obtains for Su in the limit 
of weak coupling g the usual continuum gauge field action S = ^ J d'^xF^^^. The 
staggered fermion field % corresponds to four fermion species in the continuum 
limit. The lattice model has global U(l)(g)U(l) chiral symmetry for vanishing 
fermion bare mass, mo = 0. The limit mo — 00 corresponds to the pure U(l) 
gauge theory with extended Wilson action. 

Throughout this work we consider the theory on hypercubic lattices with pe- 
riodic boundary conditions in all directions except the antiperiodic boundary 
condition of the fermion field in one ( "time" ) direction. We call such lattices 
"toroidal". 

Prom the accumulated numerical evidence we know that for any fixed value 
of mo this lattice gauge theory has a line of phase transitions in the (/?, 7)- 
plane between the strong coupling confinement phase and the weak coupling 
Coulomb phase. On finite toroidal lattices it has a two-state signal (gap Aep) 
in (cosOp) in a broad region of 7- values, including 7 = (Wilson action) 
[8,9,12,18-20]. 

These metastability phenomena are similar to those in the pure U(l) lattice 
theory. The recent studies of that theory on spherical lattices strongly suggest 
that in spite of non- vanishing Aep the order changes at 7 = 70 — 0, 70 
probably being slightly positive, and is of 2"^^ order for 7 < 70 [14,15]. In both 
theories the two-state signal on toroidal lattices weakens with decreasing 7, 
but is present until the large autocorrelation time at large negative 7 makes 
simulations prohibitively expensive. Thus on toroidal lattices the two-state 
signal on the phase transition line cannot be avoided. However, an analogy 
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Fig. 1. Expected schematie position and nature of the confinement-Coulomb phase 
transition sheet in the three-dimensional phase diagram of compact lattice QED 
with staggered fermions. The sheet of phase transitions is of 1^* (2°'^) order above 
(below) the line X^oXo. The point Xq might be at 7 = — oo. The 7 axis extends 
from large negative 71 to 72 ^ 0.5. 

with the pure U(l) theory suggests that in the full theory for any mo > the 

transition becomes 2°^^ order at some finite negative value of 7. Fig. 1 shows 
the conjectured sheet of phase transitions in the thermodynamic limit. 

For fixed 7 (and P at the transition) the two-state signal increases with de- 
creasing mo. Thus the supposed region of 2°^^ order shifts to lower values of 
7. Whether the order changes at finite 7 also at mo = (as suggested in fig 1 
by the position of the point Xq) is a challenging question. 

The continuum limit of the model can be considered at any point of the 2^^^^ 
order surface below the line X^oXo, as well as at the line itself. We expect 
that for mo > this limit will be in the universality class of the pure gauge 
theory for the following reason: In the continuum limit the lattice constant a 
is determined in physical units (e.g. by fixing some mass in physical units). 
The fermion bare mass must be expressed in such units too, 



For nonzero fixed mo the physical fermion mass diverges in the a ^ at 
nonzero mo, and fermions decouple. Thus for mo > we expect to find the 
same non-Gaussian universality class as in the pure gauge theory investigated 



mo = a mg 



phys 



(6) 



in [17]. 
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A continuum limit with finite itiq requires mo oc a — > 0. Furthermore, if 
a chiral symmetric theory itlq^^^ — should be constructed, then mo must 
approach zero faster than a. 

It should be pointed out that, if 7 < 0, the action (2) does not satisfy reflec- 
tion positivity. Of course, reflection positivity is a sufficient but not a necessary 
condition for unitarity, so its absence by itself does not invalidate the investi- 
gation of a lattice model for the purpose of quantum field theory. However, it 
makes it more demanding, as possible sources of unitarity violation must be 
identified and their scale determined. If they are on the cut-off level, they do 
no harm, at least if the cut-off can be removed. 

In our case the fermionic part of the action is reflection positive and the 
problem is due to the gauge field action at 7 < 0. Numerical investigation 
of that theory on spherical lattices [14,15] strongly suggest that the critical 
behaviour at 7 < and 7 = belongs to the same universality class. The 
latter case (Wilson action) is reflection positive, however. Thus the unitarity 
violations, if any, seem to be conflned to small (lattice) distance and do not 
appear on the physical scale. Of course, this expectation should be tested 
whenever possible. 



3 Dynamical fermions: Two-state signal 

Here we present our results of dynamical fermion simulations on 6^ toroidal 
lattice at various values of 7 and mo, concentrating on 7 < 0. The phase 
transition in the compact QED with light staggered fermions and 7 < has 
been studied last quite some time ago [9,12], to our knowledge. The two- 
state signal was observed down to large negative values of 7 = —1.0. The 
availability of the hybrid Monte Carlo fermion algorithm and larger computer 
resources suggest the reconsideration of this transition, as it may be crucial 
for understanding of QED at strong coupling. We work on toroidal lattices 
like [9,12]. 

Expecting that the lattice topology contributes substantially to the two-state 
signal, we put emphasis on the comparison of this signal with the analogous 
signal in the pure U(l) gauge theory on toroidal lattices. Therefore we study ep 
in both theories under comparable conditions . The values of Aep are obtained 
in the following way: At a chosen value of 7 we made simulations at several 
values of f3 in the vicinity of the expected phase transition. Inspecting the 
time evolution of ep the f3 value with clearest two- state signal was chosen. The 
values and errors of Aep have been obtained from a flt of the corresponding 
histogram with two Gaussians at the distance Aep. 
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Fig. 2. Time evolution of the plaquette ep at 7) = (2.08, —1.3) and mo = 0.1 on 
6^ lattice. Each point represents an average over 100 measurements. 

Table 1 

Positions of the transition and the size of the gap Aep on 6^ lattice in full QED. 



7 


mo 


/3c 


Aep 




7 


mo 


/5c 


Aep 





0.1 


0.883 


0.072(4) 




-0.2 


0.01 


0.970 


0.054(2) 


-0.2 


0.1 


1.015 


0.047(3) 




-0.2 


0.02 


0.977 


0.053(2) 


-0.5 


0.1 


1.240 


0.034(2) 




-0.2 


0.04 


0.985 


0.054(2) 


-1.0 


0.1 


1.730 


0.020(3) 




-0.2 


0.07 


1.000 


0.050(4) 


-1.3 


0.1 


2.080 


0.018(3) 




-0.2 


0.10 


1.015 


0.047(3) 


-1.5 


0.1 


2.350 


0.013(2) 











The two-state signal in full QED is clearly observable on 6^ lattice at least 
until 7 = —1.5. The values of Aep and the corresponding values of Pc on this 
lattice are collected in table 1. 

As an example we show in fig. 2 the time evolution of ep in a run at (/?, 7) = 
(2.08, —1.3) and mo = 0.1. The gap between both states is Aep ~ 0.018. The 
observed tunnelling period indicates that we need very long observation time: 
the figure represents 2.2 x 10^ hybrid Monte Carlo trajectories and the lifetime 
of the metastable states is of the order of 4 x lO'^ trajectories. This explains 
why the two-state signal was not observed for the same parameters in [12]. It 
makes also clear that increasing the lattice or further substantial decreasing 
of 7 is prohibitively expensive even today. 

The question is whether the two-state signal is a finite size effect or is caused 
by a genuine 1'^* order transition. Therefore we want to point out the small 
size of Aep: At 7 = —1.3 and mo = 0.1 it is smaller than that in the pure 
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7 mo 
Fig. 3. a) Plaquette discontinuity Aep on 6^ lattice as a function of 7 in the full 

QED for niQ = 0.1 (circles) and in the pure U(l) gauge theory (crosses). At 7 = —0.2 
also the value of Aep at rriQ = 0, obtained by extrapolation, is shown, b) Plaquette 
gap Aep on the same lattice as a function of mo in the full QED for 7 = —0.2. A 
linear extrapolation to mo = is indicated. 

U(l) gauge theory at 7 = —0.5, where the 2'^'* order behaviour is clearly seen 
on spherical lattices [14,15]. 

Similar observations can be made for mo = 0.1 in a broad interval of 7 between 
—0.2 and —1.5. As is seen in fig. 3a, there the value of Aep is smaller than 
that of the pure U(l) gauge theory at 7 = 0. We have also observed that Aep 
in the pure gauge theory roughly agrees with that in the model with fermions 
of small mass, if compared at identical /3. 

The investigation of the mo dependence of Aep becomes very expensive if one 
approaches small mo, even on the small lattice of size 6^. We have performed 
such a study at 7 = —0.2, with mo as small as mo = 0.01. The results are 
shown in fig. 3b. The increase of Aep with decreasing mo is moderate and there 
is no indication of a sudden change at very small mo- A linear extrapolation to 
mo = 0, suggested by the shape of the data, gives Aep ~ 0.055. A quadratic 
extrapolation gives an insignificantly smaller value. This value still does not 
exceed the size of the gap in the pure gauge theory at 7 = on the same 
lattice. 

The comparison of the gap values on 6^ lattice in the theory with dynamical 
fermions and in pure gauge theory would suggest that even at mo = the 
phase transition may be already of the 2"^*^ order at 7 = —0.2 (since we have 
arguments, that the pure gauge theory is really second order already at 7 = 0). 
In spite of our expectation that the gap is caused by gauge fields such a 
reasoning might be too naive. Nevertheless, it makes plausible our hypothesis 
that the phase transition at mo = becomes 2^^^^ order at some finite 7, though 
possibly smaller than 7 = —0.2. 
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4 Quenched approximation 



4-1 Simulations and fermionic observables 



The quenched calculations have been performed at 7 = —0.2 on large toroidal 
lattices, typically of the size 16^32, at some points up to 20^40. We used the 
configurations produced for the measurement of the gauge-ball masses [17] . We 
avoided the region very close to the phase transition because of the two-state 
signal there. This explains a gap in /^-values in the data and figures below. In 
this way we also avoided significant finite size effects. 

The transition point in f3 as determined from the scaling of the gauge-ball 
masses in the confinement phase at 7 = —0.2 is [3^ = 1.1607(3) [17]. Our 
following analysis is crucially dependent on a good knowledge of and will be 
based on this value. We have investigated in both phases the chiral condensate, 
(7 and TT susceptibilities and some meson masses. The names of mesons - 
fermion-antifermion bound states - are chosen in analogy with QCD. 

Following standard notation we define the chiral condensate 

(xx) = (trM-^) (7) 



and measure it with a Gaussian noise estimator. Wc also determine the loga- 
rithmic derivative Rt^ of the chiral condensate [21,22] 



^ din mo 



mo d{xx) 



I3,K 



ixx) dmo 



(8) 



I3,K 



This may be rewritten as a ratio of zero momentum a and tt meson propagators 
(susceptibilities) 



djxx) 



Ca{p = 0) , 



ixx) 

mo 



= a(p = o). 



(9) 



so that 



_ a(p ^ 0) 



The first of eqs. (9) holds in the quenched approximation, if only the connected 
part of the a meson propagator is considered. The second is a Ward Identity 
[23] . The meson operators are defined below. 
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We expect the validity of the equation of state in analogy with magnetic 
systems, [21] 

mo = {xx)^f{x), x=-^^^pj^, t = pc-p- (11) 



The suffix x of a-dded to the usual magnetic exponent /3 to avoid confusion 
with the coupling (3. 

A non-vanishing chiral condensate requires a positive zero of the scaling 
function f{x). The predicted scahng behaviour is 



ixx) = ' 



{t/xof^ for P < Pc, mo = , ^^^^ 
(mo//(0))V'5 for /3 = /3„ mo ^ . 



The equation of state (11) predicts that the ratio i?^ depends only on the 
scaling variable x, 



The timeslice operators for the mesons (projections to zero momentum) are 
given by 

0'\t) = j:s%x,,X^,t: (14) 



with the sign factors s't^, x — x,t, given in table 2. We measure their correla- 
tion functions using point sources and considering only the connected parts. 
Further details can be found e.g. in [24]. 



4-2 Chiral condensate 



The data for the chiral condensate are shown in figs. 4 and 5. In both phases 
the data turned out to be almost independent of the volume for the 16^32 and 
larger lattices we used. Figure 4 exhibits the chiral condensate as a function of 
mo for different f3 values. To extrapolate the condensate into the chiral limit 
we made for fixed (3 in the broken symmetry phase a fit with the ansatz 

(XX) (mo) = (xx)o + Amo + BrriQ In mo , (15) 
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Table 2 

List of measured mesonic operators in the lattice QED with staggered fermions. The 
continuum quantum numbers J^^ and the name of the corresponding QCD particle 
are included. The sign factors s'^ use the phase factors -q^x = (—1)^1"' 

= (_i)^M+i+ -+^-4 and Ex = (-l)^i+ - +^4^ xhe index s and a at the J^"^ are for 
the singlet and adjoint representation of the flavor symmetry group. The notation 
follows loosely [24,25]. 



i 






particle 


1 


1 


f\A — \- 


^ (fo) 


2 








3 


Vkx^xCkx 


K- 




4 


V4:x'nkx^xCkx 




h 

p{2) 



motivated by the chiral perturbation theory [26]. In the Coulomb phase we 
use the ansatz 

{xx){mo)^{xx)o + C{moy . (16) 

The value {xx)o is our estimate for the chiral condensate in the chiral limit. 

In the confinement phase we get in fact consistent results by both expressions. 
The exponent z is nearly independent of j3 , being z ^ 0.8 in the confinement 
phase, significantly different from 1. In the Coulomb phase the value is z ~ 0.7 
close to the transition and rises with /3 to z fa 0.8 at /3 = 1.2. 

The results for {xx)o ^-^e shown in fig. 5 as a function of f3. A fit ct^^ with 
P approaching the critical coupling /3c from below results in the value of the 
magnetic exponent ~ 0.33. The scaling region has been assumed to he 
between P — 1.13 and P — 1.16. Reducing this interval tends to decrease P^. 

4-3 Susceptibility ratio and equation of state 

The study of J?^ as suggested in [21] allows for a more sophisticated scaling 
analysis. As can be seen from (13), -R^r should be equal to 1/5 for P = Pc, 
as long as ttiq is small enough. In the confinement phase with broken chiral 
symmetry vanishes in the chiral limit, as can be seen from the definition of 
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Fig. 4. Chiral condensate on the 16^32-lattice as a function of mo for different /3. 
The hnes are our extrapolation into the chiral limit according to (15) (confinement 
phase) and (16) (Coulomb phase). 




Fig. 5. Chiral condensate on 16^32-lattice as a function of p. We show the data for 
different mo and our extrapolation into the chiral limit (circles). The fit has been 
done according to (12). 
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Fig. 6. i?7r as a function of mo for different (3 on 16^32 lattice. The lines connecting 
the points at the same (3 are drawn to guide the eye. 

i?7r- In the symmetric phase a and tt meson masses are degenerate, and thus 
= 1 in the chiral limit. The curve of constant i?^ separates two groups of 
curves. It provides an estimate of the critical coupling and 1/d. 

In both phases the data for i?^ are again almost independent of the volume 
for the 16^32 and larger lattices, thus we show in fig. 6 the data for 16^32 
lattice only. Since no simulation close to the critical coupling could be done, 
our estimate for 6 with this method is bound to have a large uncertainty. From 
the value of i?,r at mo = 0.05 for the /3-values closest to the phase transition 
in each phase we find S — 1.7 — 3.3. The mean-field value would he 6 — 3. 

In the symmetric phase all curves should run towards = 1. Our data do 
not show this behaviour in a clear way. However, in particular at small mo, we 
observe large statistical errors near the transition point. Also, the expectations 
of chiral symmetry restoration by mass symmetry of parity partners may not 
hold in the quenched approximation. 

To check the validity of the equation of state we also investigate the depen- 
dence of i?^ on the scaling variable x in the confinement phase, which is pre- 
dicted by (13) to be the only parameter. This is done by determining by 
the requirement that all data for R^^ lie on a single curve depending only on x. 
We get the best results for (3^ = 0.32(2) (fig. 7a), which is in good agreement 
with the result obtained from the scaling of the chiral condensate. Varying 
Pc — 1.1607(3) within its error bars neither improves the independence on x 
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nor changes (3^ significantly. 



.40 r 




Fig. 7. (a) Rt!-{x) and (b) f{x) = mo{xx)~ as functions of the scahng variable 
X = tll(xX)-^lf^^ on 16^32 lattice for the exponents (i^ = 0.32 and 5 = 1.82. The 
method to obtain these plots and the curve in (b) are discussed in the text. 



To determine the exponent 6 more precisely, we use the fact, that f{x) should 
be a function of x alone. With (11) as an indirect definition for / and taking 
our best estimate for = 0.32 we tune 6 so that mo{xx)~^ depends only 
on X. As shown in figure 7b, this leads to an excellent description of the 
data with 6 = 1.82(10). The result for 6 is substantially more precise than 
the above attempt to use i?^ directly. If this analysis is done with a slightly 
different (3^, the obtained results are much worse, and no unique value of 5 
could be obtained. The scaling function f{x) can be very well parameterized 
as f{x) = 0.207 (4.68 - x)^-^^ which is drawn in fig. 7b. 

Both exponents determined by this method, (3-^ = 0.32(2) and 5 = 1.82(10), 
differ from corresponding mean field values. The results therefore indicate that 
this model might have a non-Gaussian fixed point. They may be considered as 
a first hint that the full compact QED is nontrivial. However, we cannot give a 
reliable estimate of the systematic errors. We note that i/3^(5 + 1) = 0.23(2), 
which does not agree with the value u = 0.35 as would be expected from the 
hyperscaling relation. 



4-4 Meson spectrum 

The masses have been determined with the usual cosh-fits to the connected 
propagators [27]. Each of the four (cf. table 2) correlations functions C^t) 
contains signals from both parity states (with the only exception of the pion 
in the channel i — 2) and therefore is fitted to the form 
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Fig. 8. The figure shows vs. mo for /? = 1.10 (circles), 1.14 (crosses) and 1.158 
(squares); the fit (mjr)^ oc (mo)^ gives values of p close to one. 

n+ 

B+ + {-ifB- + (e"^"* + e-^" (^*-*)) 

n=l 

+ (-1)*^^;; (e-^"* + e-^"(^*-*)) . (17) 

n=l 

Significant mass determination could be obtained in the o", 7r*^^\ ai, p*^^^ and 
p*^^'' channels. For them, respectively, the numbers of energy levels (n^,n~) 
considered were in the confinement phase (1,1), (0,2), (1,1), (2,1), (1,2) and 
in the Coulomb phase (2,1), (0,2), (2,1), (2,1), (1,2). The lowest energy level 
was interpreted as the mass of the corresponding particle state. Statistical 
fiuctuations prevent a reliable determination of tt^^^ and hi. 

4-4- 1 Confinement phase 

In the confinement phase finite volume effects become small already at lattice 
sizes 12'^24; we therefore discuss meson masses obtained on 16^32 lattice only. 
We note that in this phase the values of both parameters B"^ obtained by the 
fits are consistent with zero. 

Close to the chiral hmit one expects 

(m^)^ oc mo . (18) 
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In quenched studies there may be corrections due to so-called chiral logarithms 
[28-30], thus we fit our results (fig. 8) to (m^)^ oc (mo)^. We find values of p 
close to one (between 1.02 and 1.13 close to the phase transition). 

The mo dependence of the non-Goldstone state masses has been parameterized 
linearly, m oc mQ + const, as usual. Fig. 9 gives the values of a and p*-^-* as they 
emerge from the extrapolation to mo = 0. The mass of p*^^^ is compatible with 
that of p^-^^. The oi is heavy (m > 1.5) and seems not to scale at the phase 
transition. 

The data is not good enough to warrant a fit to a critical exponent. However, 
assuming the non-Gaussian value u = 0.35 obtained in [14,15,17] and the 
knowledge of /3c, a fit m = ci'^ of the p^^^-mass is shown in the plot fig. 9 
to be compatible with the data. The a has mass values of the same order of 
magnitude as the p whereas the ai is definitely heavier by almost a factor of 
2. 

Comparing the p and a with the the gauge-balls (determined from the same 
runs in [17]) we find meson masses quite close to the Ti'" gauge-ball group. 
For example, the amplitude of p is c = 6.0(2) and the corresponding Ti'~- 
gauge-ball amplitude is 5.4(3). 
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Fig. 10. The tt mass in the Coulomb phase at /3 = 1.17 plotted vs. l/L^ for various 
rriQ. 

4-4-^ Coulomb phase 

In the Coulomb phase we expect and we do find sizeable finite volume depen- 
dence, presumably caused by massless photons and light fermions. The mesons 
are therefore expected to be resonances. This is like in the recent study of 
gauge-ball masses in the Coulomb phase [17]. The volume dependence there 
was explained as a signal of multi-photon states, but a reliable quantitative 
extrapolation to the thermodynamic limit was not possible. For mesons the 
situation is even worse. In particular the Bohr radius, proportional to the in- 
verse fermion mass (which we have not determined), might be much larger 
than the lattice. Therefore our investigation of the spectrum in the Coulomb 
phase has only an exploratory character, checking whether unexpected phe- 
nomena might be seen in the Coulomb phase close to the phase transition to 
the confinement phase. 

Fig. 10 indicates, that the size dependence may follow a linear 1/Ls behaviour. 
Using the correspondingly extrapolated values for the pion (para-onium) mass, 
we obtain fig. 11. For increasing fermion masses one expects itIt^ ~ 2mo, i.e. 
twice the mass of the constituents. For small constituent mass, one would 
expect oc rriQ like the positronium spectrum. However, it is not completely 
excluded that the massless QED in the Coulomb phase at strong coupling 
shows some nonperturbative phenomena due to IR singularities. 
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Fig. 11. The values of rriTr extrapolated to infinite volume against mo. We also show 
the results of fits (19) with (solid line) and without constraint A = 0. 

Therefore we have fitted the mass values with a simple three-parameter curve 

A + Brno + 2ml 



C + mo 



(19) 



either keeping A fixed to (solid line in fig. 11) or letting A free (dashed line, 
A = 0.1(1)). We found, that the data are compatible with both, vanishing 
and non- vanishing values of 171.^ at mo = 0, and the conventional expectation 
of vanishing meson masses is in agreement with the data. The non-vanishing 
value would imply a decoupling of the corresponding state in the continuum 
theory if we require that a continuum limit should be taken at /? ^ /?c and 
mo —>■ 0. 

Comparing the mass values for all states as obtained on a fixed lattice size 
16^32 we observe an approximate mass degeneracy of the chiral partners a 
and TT even at finite mo, consistent with chiral symmetry in the mo = limit. 
This is in agreement with the discussion of i?^ approaching 1. Also p^"^^ and 
oi masses are close to each other and to the a and tt masses. Thus the data 
do not show any significant fine or hyperfine splitting. 

This explorative study of the spectrum indicates that much better data and 
understanding of finite size effects are required to obtain results for the spec- 
trum in the Coulomb phase which would bring further insight. 
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5 Summciry 



We have studied the compact U(l) gauge theory with staggered fermions on 
toroidal lattices on both sides of the phase transition line in the (/?, 7)-plane. 



• With dynamical fermions the observed gap Aep on 6'^-lattices increases 
with respect to the pure gauge theory, decreases with decreasing 7 but does 
not vanish in the considered domain —1.5 < 7 < 0. At 7 = —0.2 and for 

masslcss fermions it is of the same size as for the pure gauge theory at 
7 = 0. Since there the gap disappears on sphere-hke lattices wc conjecture, 
that also the fermionic theory may in fact have a critical phase transition 
at negative 7. 

• On large lattices we study the chiral condensate and the meson spectrum in 
the quenched approximation. Quenched results should not be over-interpret- 
ed. We hope, however, that these results stimulate further studies of the 
compact lattice QED with dynamical fermions. 

• The chiral condensate is determined with high precision. The main result 
of our study is the observation that the chiral condensate can be very well 
described as a function of both parameters, f3 and mo, by means of an 
equation of state with the relevant critical indices that have non-Gaussian 
values: (3 ~ 0.32 and 6 ~ 1.8. 

• In the confinement phase we observe chiral symmetry breaking and Goldsto- 
ne-boson behaviour of the pion, as expected. The mass squared obeys a 
PCAC behaviour. All masses (extrapolated to vanishing bare fermion mass) 
have a behaviour compatible with a non-Gaussian scaling at the phase tran- 
sition in (3. 

• The results for masses in the confinement phase are generally of better 
statistical quality than the Coulomb phase results. In the Coulomb phase 
the masses of the pion and other states are of comparable size. 



Wc obtained a better insight into some aspects of compact lattice QED, in 
particular concerning chiral symmetry breaking and the meson spectrum in 
the confinement phase. However, higher statistics will be necessary to improve 
the quality of some of our results in this phase and to determine the scaling 
behaviour of the meson masses. In the Coulomb phase not only considerably 
higher statistics, but also understanding of finite size effects of various origin 
is required. Here the problems are similar to those of the noncompact QED, 
except that (3c can be determined from the other, technically simpler phase. 
This might be an advantage of the compact formulation. 
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